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Proton migration is a ubiquitous process in chemical reactions related to biology, combustion, and 
catalysis. Thus, the ability to control the movement of nuclei with tailored light, within a hydro¬ 
carbon molecule holds promise for far-reaching applications. Here, we demonstrate the steering of 
hydrogen migration in simple hydrocarbons, namely acetylene and ahene, using waveform-controlled, 
few-cycle laser pulses. The rearrangement dynamics are monitored using coincident 3D momentum 
imaging spectroscopy, and described with a quantum-dynamical model. Our observations reveal 
that the underlying control mechanism is due to the manipulation of the phases in a vibrational 
wavepacket by the intense off-resonant laser held. 

PACS numbers: 33.80.Rv, 82.37.Vb, 82.50.Nd 

The rearrangement of hydrocarbon bonds via the mi¬ 
gration of a hydrogen atom can result in major de¬ 
formations of molecular architecture, and thus alter 
the molecule‘s chemical properties. Examples include 
keto-enol tautomerism where the migration of a proton 
changes an aldehyde into an alcohol. Isomerization re¬ 
actions of that kind have been the subject of numerous 
studies mn]. Of particular interest was to determine 
the so-called isomerization time, which has been mea¬ 
sured to be within several tens of femtoseconds in small 
hydrocarbons nn. The phenomenon has also been ob¬ 
served in larger molecules, such as protonated triglycine 
[7]. Tracing of the hydrogen migration from different lo¬ 
cations within the molecule has been made possible via 
isotope labeling, see e.g. [HHin]- The ability to exert 
control over the migration could lead to advancement in 
topics such as the efficiency of catalytic reactions m and 
combustion reactions regarding fuel and energy research 
na. Furthermore, light-induced control of hydrogen mi¬ 
gration may open new reaction pathways which cannot 
materialize by other means. 

Despite its direct relevance to applied chemistry, stud¬ 
ies regarding the control of the hydrogen migration pro¬ 
cess have been scarce, and have been limited to theory 
[MlIIl] for a long time. However, recently progress has 
been made in coherently controlling isomerization reac¬ 
tions using fundamental parameters of ultrafast strong- 
held laser sources. Xie, et al. varied the pulse dura¬ 
tion and intensity to explore the isomerization of ethy¬ 
lene US], and reported control of the total fragmentation 
yields of various hydrocarbons m- 

Here, we demonstrate control over the direction of the 
hydrogen migration using the electric-held waveform of 


intense few-cycle laser pulses. This approach goes be¬ 
yond earlier work on tolluene m and methanol CHI us¬ 
ing two-color pulses with a duration of tens of fs. In 
contrast, the duration of our few-cycle laser pulses is sig- 
nihcantly shorter than the time scale of the isomerization 
dynamics, therefore avoiding charge-resonance-enhanced 
ionization m occuring at large internuclear distances 
m- Moreover, the an inhuence of electron localization- 
assisted enhanced ionization on the dissociation reac¬ 
tions, recently demonstrated for acetylene theoretically 
[2T1 122] and experimentally [23] , is also avoided. In com¬ 
parison to earlier experiments using longer pulses, our 
experiment therefore allows the study of the pure prede¬ 
termination of a molecular reaction, in the sense that the 
laser held initiates a reaction but is not present during 
the actual reaction. 

The electric held of a few-cycle laser pulse can be de¬ 
scribed in the time domain as E{t) = Eo(t) cos {ut -h 0), 
where Eq (t) is the envelope, uj the carrier frequency, and 
(f is the phase between the carrier and the envelope (car¬ 
rier envelope phase; CEP). The CEP determines the elec¬ 
tric held waveform on a sub-femtosecond timescale and 
has proven to be a powerful parameter for controlling 
electron motion in various systems pni It has 

also recently been shown to apply to the control of nu¬ 
clear wavepackets m- 

In the present work, we show that the CEP of a 
few-cycle laser pulse further permits to control complex 
photochemical reactions involving structural rearrange¬ 
ments. Using a reaction microscope (REMI) [28], we ob¬ 
serve asymmetric emission of carbon ions from acetylene 
and trihydrogen ions from allene as a function of the 
CEP of the driving laser pulses, demonstrating the con- 
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trol of proton migration from one side of a molecule to 
the other side. The results are interpreted in terms of a 
quantum mechanical model where the direction of hydro¬ 
gen migration results from the phases of a superposition 
of vibrational modes manipulated by the laser CEP. 

The experimental setup will be described in detail else¬ 
where. Briely, hydrogen migration in acetylene and al- 
lene is induced and controlled by intense, CEP-stable 
4-fs laser pulses at a carrier wavelength of 750 nm and 
10 kHz repetition rate. The near-single-cycle pulses are 
obtained from the SMILE laser system at LMU, which 
is based on the Eemtolasers Eemtopower HR CEP4 am¬ 
plified laser system. The CEP of the laser is stabilized 
using the feed-forward technique [29] and measured us¬ 
ing either a stereo-ATI phase meter m (in the acetylene 
experiments) or an f-2f interferometer [31] (in the allene 
experiments). 

In the REMI, the laser pulses are focused (/ = 17.5 cm) 
into a cold gas jet of neutral hydrocarbon molecules. 
A homogeneous electric field is used to direct all ions 
generated in the laser focus onto a time- and position- 
sensitive multi-hit capable detector, providing the three- 
dimensional (3D) momentum distributions of multiple 
ions in coincidence. 

To study the strong-held control of hydrogen-bond re¬ 
arrangements in acetylene (HCCH), we focus on its iso¬ 
merization to vinylidene (CCH 2 ) via proton migration 
from one side of the molecule to the other. Vinylidene 
formation is clearly identihed by the dissociation of the 
dication into C+ and CHj fragments. The distinct sep¬ 
aration between the fragments from the acetylene and 
vinylidene breakups, which have similar mass and mo¬ 
menta, is challenging and feasible so far only by coinci¬ 
dence detection methods [28l [32] . 

Eig.[^ shows the experimental results for the vinylidene 
breakup after isomerization of acetylene. The momen¬ 
tum distribution of ions (Eig.j^a)) is nearly isotropic. 
This can be understood as a consequence of the hydro¬ 
gen migration preceding dissociation of vinylidene along 
the C-C bond [32]. Eor C+ momenta in the indicated 
range —42a.u. < Pz < —28a.u., the C+ ions arrive at 
the ion detector almost simultaneously with the coin¬ 
cident CHJ fragments, which impedes their coincident 
detection. Hence only events with momenta outside the 
range 28a.u. < \pz\ < 42a.u. are considered for evalu¬ 
ating the CEP-dependence of the C+ emission direction. 
The emission direction of the ion following the cleav¬ 
age of the C-C bond is directly indicative of whether the 
left proton in acetylene has migrated to the right side 
(along the laser polarization), or vice versa. The yields 
of C+ ions emitted to the right, R( 0 ), or left, 1 /( 0 ), as 
a function of the CEP are therefore a measure of the 
CEP control of the hydrogen migration. To quantify the 
control, we introduce the asymmetry parameter as 

A{cp) = {L{4>) - nm / m) +R{cf>)). ( 1 ) 



FIG. 1. Strong-field control of hydrogen migration in acety¬ 
lene. (a) Measured CEP-averaged momentum distribution 
for C+ ions emitted from vinylidene dications following hy¬ 
drogen migration in acetylene, induced by 4-fs laser pulses 
with an intensity I ~ 1.3 x lO^^W/cm^. The shaded area is 
cut out for technical reasons (see text), (b) CEP-dependent 
asymmetry parameter for the C^ yield shown in (a), eval¬ 
uated for a cone of 80° around the polarization axis of the 
laser. Since the absolute value of the CEP is not known from 
the experiment, the phase is given relative to the asymme¬ 
try parameter recorded for C 2 HJ recoil ions (see text). The 
solid line depicts the calculated asymmetry curve for an in¬ 
tensity of 1.5 X lO^^W/cm^, shifted along the CEP axis for 
best agreement with the experiment. 


The CEP-dependence of the asymmetry parameter 
Ac+(0), recorded for C+ ions and plotted in Eig. [^b), 
demonstrates the steering of the direction of proton mi¬ 
gration in acetylene. 

Eor allene (H 2 C 3 H 2 ) molecules, the hydrogen motion 
is monitored through the formation of trihydrogen ions 
(H^). This reaction channel is identified by the dissocia¬ 
tion of the allene dication into H^ and C 3 H+ fragments. 
Even though allene’s isomer propyne (HC 3 H 3 ) is stable 
in the neutral state and contained in low concentrations 
in commercially available allene gas bottles, H^ forma¬ 
tion has been found to result predominantly from the 
isomerization of allene [ 8 ]. 

The results for the H^ formation from allene are dis¬ 
played in Ei g. [^ The momentum distribution recorded 
for H+ (Eig. m ^a)) is nearly isotropic and therefore cor¬ 
roborates the finding that H^ is predominantly formed 
following hydrogen migration. Thus, the CEP depen¬ 
dence of the asymmetry parameter A^+ ( 0 ) demonstrates 
the steering of the direction of proton migration in al¬ 
lene. Interestingly, we find in our experiments that the 
asymmetry in the H^ formation from allene exhibits a 
pronounced sensitivity to the laser intensity [33] . 

The CEP axes in Eigs. Bb) and[^b) are given relative 
to the phase of the asymmetry parameters (see eq. 
measured for the recoil cations of either molecule, pro¬ 
duced by single ionization. This makes a comparison of 
the phases of the asymmetry parameters recorded for hy¬ 
drogen migration in acetylene and allene possible, indi¬ 
cating that Ac+ ( 0 ) and ( 0 ) assume maxima for sim¬ 
ilar but not identical CEP values. The experimental data 
is further compared to predictions from our quantum me¬ 
chanical calculations. Because the absolute CEP in the 
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FIG. 2. Strong-field control of hydrogen migration in al- 
lene. CEP-averaged Momentum distribution (a) and CEP- 
dependent asymmetry parameter (b) for HJ ions emitted 
from propyne dications following hydrogen migration in al- 
lene. The CEP-dependence of the is evaluated for a cone 
of 80° around the polarization axis of the laser and at a peak 
intensity of ~ 2.9 x lO^^W/cm^. The phase is given relative to 
the asymmetry parameter recorded for CsHJ recoil ions (see 
text). The measured asymmetry curve is compared to the one 
calculated for an intensity of 3.0 x lO^^W/cm^, shifted along 
the CEP axis by the same value as the theory curve in Eig. 

0b)- 


experiment is not known, the predicted curve for 
emission from acetylene was shifted along the CEP for 
best agreement with the experiment. The prediction for 
emission from allene was shifted by the same amount 
and lines up with the measured asymmetry. Moreover, 
the amplitudes of the calculated asymmetries agree rea¬ 
sonably well with the measured ones, lending support to 
the model described in the following. 

In general any nuclear motion, can be expressed in the 
basis of all normal modes. We use normal modes to de¬ 
scribe the molecule laser interaction and reactive coordi¬ 
nates for the description of the subsequent hydrogen mi¬ 
gration. The subset of normal modes relevant for the ini¬ 
tiation of the reaction is illustrated in Fig. For Acety¬ 
lene (Fig. [^a)) the IR-active cis-bending mode |n0) and 
the IR-inactive trans-bending mode |0m) form a 2D basis 
|nm), where m,n are the number of vibrational quanta. 

Here, the respective time evolution factor exp 
is implicitly included and the relative phase of the eigen¬ 
functions is set to compensate any phase offset relative 
to the CEP. For allene Fig. [^b) three normal modes 
(IR-active rocking and anti-symmetric bending and IR- 
inactive symmetric bending) are required to describe the 
initiation of isomerization by the laser field. 

The interaction of the molecules with the external 
laser field is calculated by solving the time-dependent 
Schrodinger equation (TDSE) using the Chebychev prop¬ 
agator. The TDSE for the neutral molecule under the 
influence of the light field reads: 

(^n + (t)) (i) = (2) 

where is the Hamiltonian of is the as¬ 

sociated dipole moment, and (t) is the nuclear wave 
function. The light-field e(t) is included in the dipole ap- 
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FIG. 3. Reactive coordinates and normal modes used in the 
theoretical description of the hydrogen migration reactions in 
acetylene and allene. The wavepacket formation in acetylene 
(allene) is calculated on two- (three-) dimensional potential 
energy surfaces (PES) along the normal modes shown in pan¬ 
els (a)((c)). The reactive coordinates (^ 1 ,^ 2 )- (^i,<f 2 ) repre¬ 
sent the CCH bond angles (shown in cyan) and are used to 
describe the wavepacket propagation leading to isomerization 
by hydrogen migration from the left to the right, or from the 
right to the left, (b) The PES of the A^H state of the acetylene 
dication, on which isomerization occurs. The binding energy 
is plotted as a function of the two reactive coordinates. The 
initial form of the molecule and the two isomerized configura¬ 
tions are indicated, (d) Same as (b) for the B^H state of the 
allene dication, relevant for the isomerization reaction which 
leads to HJ formation. 


proximation and is characterized by a full width at half 
maximum (FWHM) of the intensity envelope of 4fs. 

The excitation/ionization mechanism is explained ex¬ 
emplar ily for acetylene. The laser only addresses the 
IR-active modes |n0), in which it creates a wavepacket 
I 72 O). This wavepacket forms temporarily even if the 
laser pulse is non-resonant to the vibration, and the phase 
(\) of the vibrational wavepacket matches the laser CEP. 
For our simulations, we consider that the molecule is ion¬ 
ized at the peak intensity of the laser pulse. Without the 
ionization, the molecule would remain in the vibrational 
ground state after the laser pulse has passed. 

As the eigenfunctions of the molecular cation are 
slightly different from the ones of the neutral molecule, 
the ionization leads to population of several cationic 
modes, including the IR-inactive ones |0m). This pop¬ 
ulation process by ionization is independent of the CEP. 
The CEP-dependent part of the full wavepacket is there- 
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fore predominantly given by the basic superposition of 
the fundamental IR-active and the IR-inactive modes, 
^basic = 101) +exp“*'^ 110). For interpretation of the con¬ 
trol mechanism, it is sufficient to concentrate on 4^basic5 
because it is responsible for the asymmetric isomeriza¬ 
tion. In a second ionization step, the (full) wavepacket 
is projected onto the reactive, excited dicationic state, 
which is the lowest excited state we find to support iso¬ 
merization. 

In the second simulation step, the isomerization dy¬ 
namics in the reactive state are described on two dimen¬ 
sional PES along two reactive coordinates illustrated in 
Fig-S These coordinates describe the CCH bond angles 
leading to isomerization. To calculate the isomerization 
we project the wavepacket Tmodes? expressed in normal 
modes, onto the wavepacket V^reac? expressed in reactive 
coordinates. The projection is performed via the respec¬ 
tive eigenfunctions of the dication ground state 

V^reac ~ ^ ^ (^modes | ^i) | 4^j) \4^j) •> (3) 

hj 

where are the normal mode eigenfunctions and (pj are 
the reactive coordinate eigenfunctions. 

The isomerization is simulated through evolution of the 
prepared wavepacket on the reactive states for 480 fs with 
time intervals of 0.24 fs. The reactive states are the A^Ii 
state of the acetylene dication (Fig. [^b)) and the 
state of the allene dication (Fig. |^d)). On both reactive 
states the Franck-Condon region is not the global mini¬ 
mum. Isomerization occurs if the small barriers are over¬ 
come by the kinetic energy of the generated wavepackets 
and the isomerization minima indicated in Figsj^b) and 
(d), are reached. From these isomerization minima, dis¬ 
sociation can take place. The dissociation yield L (R) 
corresponding to the left (right) isomerization is quanti¬ 
fied by the part of the wavepacket at the corresponding 
isomerization minimum after integration over the propa¬ 
gation time. The asymmetry parameter A is calculated 
from L and R according to Eq. 

For the isomerization of acetylene to vinylidene, 
the control mechanism is illustrated by the calculated 
wavepacket propagation displayed in Fig. The 

wavepacket, that is prepared in the neutral and cation, is 
projected onto the reactive A^U state in the dication and 
propagated along the two reactive coordinates introduced 
above. The wavepacket possesses a CEP-dependent pref¬ 
erential propagation direction, contained in the basic 
wavepacket Tbasic- The propagation direction results in 
population of either one of the two indicated potential 
minima corresponding to the left or the right vinylidene 
configurations. This leads to the experimentally observed 
bias in the C+ emission direction after dissociation of the 
C-C bond, tunable with the CEP. 

The isomerization of allene to propyne occurs on the 
PES plotted in Fig. [^d) in a similar manner as the iso- 
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FIG. 4. Schematic of the control mechanism for hydro¬ 
gen migration in acetylene. Vibrational wavepackets of the 
Tbasic(O) = |0I) + 110) (top row) and the Tbasic(7r) = |0I) — 

110) (bottom row) superposition (see text), are propagated 
on the PES of the first excited state of the acetylene dica¬ 
tion along two reactive coordinates. The contour lines dis¬ 
play snapshots of the calculated nuclear wavepacket which 
leads to isomerization of the acetylene dication. The CEP of 
the few-cycle pulse determines the sign in the superposition 
and thereby controls whether the left proton moves to the 
right (top row) or the right proton to the left (bottom row). 


merization of acetylene but is more complex. The re¬ 
active coordinates for the allene isomerization are a lin¬ 
ear combination of eight normal modes. Three of these 
modes (one IR-inactive, two IR-active ones, see Fig. ic)) 
are populated in the preparation step, giving rise to a 
CEP-dependent wavepacket T = |00m)+exp“*'^(|ni00) + 
| 0 n 20 )), which allows for steering of the direction of hy¬ 
drogen migration. 

In conclusion, we have demonstrated experimentally 
that few-cycle laser pulses can be exploited to control 
the preferential direction of hydrogen migration in small 
hydrocarbon molecules. The hydrogen migration pro¬ 
cesses leading to the isomerization of acetylene to vinyli- 
dene and to formation in allene can be understood in 
terms of a nuclear control mechanism, similar to the one 
proposed in Ref. m- According to this mechanism, a 
field-sensitive nuclear wavepacket is prepared during the 
few-femtosecond interaction of the molecules with the 
laser field which predetermines the outcome of an iso¬ 
merization reaction taking place on a significantly longer 
time scale. The efficiency of the control could be im¬ 
proved by resonant excitation of the vibrational modes 
and driving ionization by even shorter laser pulses. Our 
work paves the way towards strong-held coherent control 
of the structural and functional properties of hydrocar¬ 
bon molecules. 
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